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Abstract 

> 

' We study the 4-dimensional pressureless Navier-Stokes-Poisson equations with density- 

dependent viscosity. The analytical solutions with arbitrary time blowup, in radial symmetry, 
are constructed in this paper. 



1 Introduction 



X 



The evolution of a self-gravitating fluid can be formulated by the Navier-Stokes-Poisson equations 



of the following form: 

C3 



Pt+V • (pu) =0, 

(pu) t +V • (pu ® u) + VP =-pV$ + vis(p, u), (1) 
A$(t,x) = a(N)p, 

where a(N) is a constant related to the unit ball in R N : a(l) = 2; a(2) = 2ir and For N > 3, 

a(N) = N(N - 2)V(N) = N(N - 2) — — - , (2) 

where V(N) is the volume of the unit ball in R N and T is a Gamma function. And as usual, 
p = p(t,x) and u = u(t,x) G R w are the density, the velocity respectively. P = P(p) is the 
pressure. 
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In the above system, the self-gravitational potential field $ = <&(<, x) is determined by the 
density p through the Poisson equation. 
And vis(p,u) is the viscosity function: 

vis(p,u) = v(mO) V ,u )- (3) 
Here we under a common assumption for: 

Kp) = (4) 

and k and 9 > are the constants. In particular, when = 0, it returns the expression for the u 
dependent only viscosity function: 

vis(p, u) = kAu. (5) 

And the vector Laplacian in u{t, r) can be expressed: 

N-l N-l 
An = u rr H u r — u. (6) 

The equations ([TJ) i and (TT])2 (vis(p,u) ^ 0) are the compressible Navier-Stokes equations with 
forcing term. The equation (JTJ)3 is the Poisson equation through which the gravitational potential 
is determined by the density distribution of the density itself. Thus, we call the system ([TJ the 
Navier-Stokes-Poisson equations. 

Here, if the vis(p,u) = 0, the system is called the Euler-Poisson equations. In this case, the 
equations can be viewed as a prefect gas model. For N = 3, ([T]) is a classical (nonrelativistic) 
description of a galaxy, in astrophysics. See [2], [7] for a detail about the system. 

P = P(p) is the pressure. The 7-law can be applied on the pressure P(p), i.e. 

P{p) = Kp<±£-, (7) 

7 

which is a commonly the hypothesis. The constant 7 = cp/c„ > 1, where cp, c v are the specific 
heats per unit mass under constant pressure and constant volume respectively, is the ratio of the 
specific heats, that is, the adiabatic exponent in |7|l. In particular, the fluid is called isothermal if 
7 = 1. With K — 0, we call the system is pressureless. 

For the 3-dimensional case, we are interested in the hydrostatic equilibrium specified by u = 0. 
According to [2], the ratio between the core density p(0) and the mean density p for 6/5 < 7 < 2 
is given by 

f ~ 3 -ii*)) (8) 



p(0) v z 



z=z a 



where y is the solution of the Lane-Emden equation with n = 1/(7 — 1), 

y(z) + -y(z) + y(z) n = 0, y(0) = a > 0, y(Q) = 0, n = (9) 
z 7 — t 



Analytical Blowup Solutions 



3 



and zq is the first zero of y(zo) = 0. We can solve the Lane-Emden equation analytically for 



Vanal (■2-) 



l-±z 2 , n = 0; 
sin z 

j n = I; 

z 

1 , 

n = 5, 



(10) 



I v/l + 2 2 /3 

and for the other values, only numerical values can be obtained. It can be shown that for n < 5, 
the radius of polytropic models is finite; for n > 5, the radius is infinite. 

Gambin [4] and Bezard |ij obtained the existence results about the explicitly stationary solution 
(u = 0) for 7 = 6/5 in Euler-Poisson equations: 

^{^)"\i + Avr'\ (id 

where A is constant. 

The Poisson equation JT])3 can be solved as 



*(t,x)= G(x-y)p(t,y)dy, (12) 

JR N 

where G is the Green's function for the Poisson equation in the TV-dimensional spaces defined by 



G{x) = < 



\x\, N = l; 

\og\x\, N = 2; (13) 



In the following, we always seek solutions in radial symmetry. Thus, the Poisson equation {T])3 is 
transformed to 

t-N-^rrfax) + (N -l)r N - 2 $ r =a(N)pr N -\ (14) 

*r = 2t TTT r p(t,s)s N - 1 ds. (15) 
r Jo 

Definition 1 (Blowup) We say a solution blows up if one of the following conditions is satisfied: 

(1) The solution becomes infinitely large at some point x and some finite time To; 

(2) The derivative of the solution becomes infinitely large at some point x and some finite time To. 

In this paper, we concern about blowup solutions for the 4-dimensional pressureless Navier- 
Stokes-Poisson equations with the density-dependent viscosity. And our aim is to construct a 
family of such blowup solutions. 

Historically in astrophysics, Goldreich and Weber [5] constructed the analytical blowup solution 
(collapsing) of the 3-dimensional Euler-Poisson equations for 7 = 4/3 for the non-rotating gas 
spheres. After that, Makino [7] obtained the rigorously mathematical proof of the existence of 
such kind of blowup solutions. And in [3j, we find the extension of the above blowup solutions to 
the case . In [S], the solutions with a from is rewritten as 
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For N > 3 and 7 = (2N - 2)/N 



y(^) N/iN - 2 \forr<a(t)Z„ &{t) 
0, for a(t)Z^ < r. a[ > 



a(t)-- 



X 



a(tW 



— , a(0) = a > 0, d(Q)= 01, 



(16) 



i)(z)+^j^y(z)+ {2 ^ N) 2)K y( z f /{N ~ 2) = M, 2/(0) = a > 0, v(0) - 0, 

where /x = [iV(iV - 2)A]/(2JV - 2)Jf and the finite Z^ is the first zero of y(z); 
For N — 2 and 7=1, 

p( t , r ) = JL e «W«M), u (t, r ) =^L 



a(t) 



a(i) a 

o(t)= T-r- , a(0) = a > 0, d(0)= a x ; 

a(f) 

1 27T ^ 

j/(z)+-y(z)+-^ e = n, ?/(0) = a, y(0) = 0, 

2 it 



(17) 



where if > 0, fi = 2X/K with a sufficiently small A and a are constants. 

For the construction of special analytical solutions to the Navier-Stokes equations in R N or 
the Navier-Stokes-Poisson equations in i? 3 without pressure with 6 = 1, readers may refer Yuen's 
recent results in [TU], [H] respectively. 

In this article, the analytical blowup solutions are constructed in the pressureless Euler-Poisson 
equations with density-dependent viscosity in R with 6 = 5/4 in radial symmetry: 



p t + Up r + pu r + -pu = 0, 

r 

o?(4)p f r 3 3 

p(u t + uu r ) = 7 — / p(t, s)s 3 ds + [np 5 / 4 ] r u r + (np b / A )(u rr + -u r ^u), 
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(18) 



in the form of the following theorem. 



Theorem 2 For the A- dimensional pressureless Navier-Stokes-Poisson equations with 9 
radial symmetry, H18\) , there exists a family of blowup solutions, 

1 ....... -G 



{ 



P(t,r) 



T-Ct ' 



y(z)+-y(z) + ^y(z) = 0, 2/(0) = a, y(0) = 0, 
z kC 



5/4, in 



(19) 



where T > 0, k > 0, C > and a are constants. 
And the solutions blow up in the finite time T/C . 



2 Separable Blowup Solutions 

Before presenting the proof of Theorem[2] we prepare some lemmas. First, we obtain the solutions 
for the continuity equation of mass in radial symmetry (|18| 1 . 
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Lemma 3 For the 4~ dimensional conservation of mass in radial symmetry 



p t + up r + pu r + -pu = 0, (20) 



there exist solutions, 



where T and C are positive constants. 



Proof. We just plug (21]) into ([20]). Then 

3 

p t + p r u + pu r + -pu 

J-4)(-C)j/(^t) 4 MT±ci) 3 y(T±ci)r(-l)(-C) 
(T-Ct) 4 (T-Ct) 4 (T-Ct) 2 

MT^ct) 3 y(T^ct) 1 (-C) v(t^ci)' (-C) 3 y(T^ct) 4 (-C) 
(T-Ct) 4 T-CtT-Ct (T-Ct) 4 T-Ct r (T-Ct) 3 T-Ct 

(T-Ct) 4 (T-Ct) 6 r 

Wy(T^r t ?y(r/(T - Ct)) Q/(^) 4 SC^^) 4 
(T-Ct) 6 T (T-Ct) 4 (T-Ct) 4 

= 0. 

The proof is completed. ■ 

Besides, we need the lemma for stating the property of the function y(z). The similar lemma 
was already given in Lemmas 2.1, [3], by the fixed point theorem. The proof is similar, the proof 
may be skipped here. 



Lemma 4 For the ordinary differential equation, 



m + - z y(*) + f^) 4 = °' (22) 

2/(0) = a > 0, 2/(0) = 0, 



where a (4), C and k are positive constants, 

has a solution y(z) £ C 2 provided that y(z) C [a, 0]. 



Here we are already to give the proof of Theorem (2] 
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Proof of Theorem 2. From Lemma|31 it is clear for that (fTO)) satisfy (fT5|) i . For the momentum 
equation (fT5|) 9. we get, 



p(u t + uu r ) H — [ p{pi s)s 3 ds — [np b ^\ r u r — Kp 5 ^ 4 (u rr H — u r ^u) 



, (-co „ (-c) 



(T-Ct) 2 T-Ct T-Ct 



a(4)p 



s \4 



'.T-Ct' 3 



r 3 7 (T-Ct) 4 ' 





(23) 
(24) 



1 / r 

K {T-Ct) iV \T-Ct 

a(4)p 7 y(T^t) 4 



-i 5/4 



T-Ct 



r 3 J {T-Ct) A 
o 



S 3 C?S ft 



1 



(T-C*f) 4W \T-Ct 



1/4 



1 (-C) 



(T - Ci) 4 T-Ct T-Ct 



a(4)p fy(T^ct) 4 



r 3 y {T-Ct)' 
o 

(4)p ?2/(t^ct) 4 



T-Ct (T - Ci) 4 (T - C*i) 2 



r 3 y (T-Ct)' 
o 



s 3 ds 



5Cftp ■ / r 



(T - Ci) 3 
P 

(T - Ci) 3 
P 



5Cfty( 



5Cfty( 



{T-Ctf V \T-Ct 



T-Ct' r 3 (T - Ct) J T-Ct 
o 

r/(T-Ct) 

a(4) /" „( B )Vd* 



T i T l ct ) 3 



Q 



(T-Ct) 3 ^ \T-Ct 
And denote 



q{ t~c1 ] = Q{z) = 5C >^)+^r 1 </<*> : 



(25) 
(26) 

(27) 

(28) 

(29) 
(30) 



(31) 



Differentiate Q(z) with respect to z, 



Q(z) 

= hC Ky(z)+a(A)y(z) 
3 



4 1 *=^/y(.)Vd. 



5Cfty(z 

z z z 



—Q(z), 

z 



(32) 
(33) 

(34) 
(35) 
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where the above result is due to the fact that we choose the following ordinary differential equation, 

y{z) + -z V{z) + 5C^ V{Z) =°- (36) 
y(0) = a > 0, y(Q)= 0. 

With Q(0) = 0, this implies that Q(z) = 0. Thus, the momentum equation (fTg|) ? is satisfied. 
Now we are able to show that the family of the solutions blow up, in the finite time T/C. This 
completes the proof. ■ 

The statement about the blowup rate will be immediately followed: 

Corollary 5 The blowup rate of the solution M9\) is 

lim p(t,0)(T -Ctf >0(\). (37) 

t-*T/C- 

Given that the sign of the constant C in (fTT|) is changed to be negative, the below corollary is 
clearly shown. 

Corollary 6 For the A- dimensional pressureless Navier-Stokes-Poisson equations, with 8 = 5/4, 
in radial symmetry, 118\) . there exists a family of solutions, 

p ^ r) = jr^w y (^y> u ^ r) = Y^ct r; m 

j/(z)+-y(z)+^y(z) 4 = 0, y(0) = a > 0, y(0) = 0, 
where T > 0, k > 0, C < and a are constants. 

Remark 7 Besides, if we consider the 4- dimensional Navier-Stokes equations with the repulsive 
force in radial symmetry, 



3 

p t + Up r + pu r + -pu = 0, 



r 



«(4V rr„n „^3,„ , r.^5/4u. , k/4v„. , 3 3 i:]<)] 



p (u t + uu r ) = H 3 — /J* p(t, s)s 3 ds + [np 5 / 4 ]u r + (np b ^)(u rr + -u r ^u), 

the special solutions are: 



P( f > r ) = /jr-W ^(ra) ,u(t,r)=- — -/ 
y(z)+-m ~ l§V (^) 4 = 0, 2/(0) = a, y(0) = 0. 



-C 

\T-CIV"\> ! "' ; ' ' T-Cf (4Q) 
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